Roton immiscibility in a two-component dipolar Bose gas 
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We characterize the immiscibility-miscibility transition (IMT) of a two-component Bose-Einstein 
condensate (BEC) with dipole-dipole interactions. In particular, we consider the quasi-two di- 
mensional geometry, where a strong trapping potential admits only zero-point motion in the trap 
direction, while the atoms are more free to move in the transverse directions. We employ the Bo- 
goliubov treatment of the two-component system to identify both the well-known long-wavelength 
IMT in addition to a roton-like IMT, where the transition occurs at finite-wave number and is 
reminiscent of the roton softening in the single component dipolar BEC. Additionally, we verify the 
existence of the roton IMT in the fully trapped, finite systems by direct numerical simulation of the 
two-component coupled non-local Gross-Pitaevskii equations. 

PACS numbers: 



I. INTRODUCTION 

The phenomenon of Bose-Einstein condensation is 
characterized by the presence of long-range phase co- 
herence. Interestingly, phase coherence can persist in a 
two-component Bose-Einstein condensate (BEC), result- 
ing in an overlapping or miscible mixture of distinguish- 
able components. This miscible state, however, is only 
stable for a certain range of interaction parameters and 
densities. Indeed, the system can be driven to immisci- 
bility by modulating the density of the two-component 
BEC, by tuning the relative strengths of the interspecies 
and intraspecies interactions or by altering the geometry 
of the trap in which the system is held [I]-|6]. 

To date, most of the experimental work on such sys- 
tems has been performed using alkali atoms that interact 
predominantly via short-range potentials |7-9 . For such 
interactions, the transition to immiscibility is character- 
ized by the parameter A = 511322/512 ~ 1) where gn and 
322 are the intraspecies interaction couplings and 1712 is 
the interspecies interaction coupling. Having A > im- 
plies the stability of the miscible state and A < implies 
an unstable miscible state. Thus, the transition to im- 
miscibility can be seen to originate from the competing 
strengths of the interspecies and intraspecies contact in- 
teractions. The immiscible-miscible transition (IMT) is 
characterized by A = 0, where the interspecies and in- 
traspecies interactions are balanced [lOl QT] . The pres- 
ence of an external trapping potential, however, relaxes 
this criterion as the trap introduces an additional energy 
cost for the components to spatially separate . 

Recently, much effort is being directed towards creat- 
ing quantum degenerate gases of atoms and molecules 
that possess non-negligible dipole moments and thus 
can interact via both short-range and dipole-dipole in- 
teractions, which are long-range (oc 1/r 3 , where r is 
the distance between the dipoles) and anisotropic in na- 
ture. Already, experimentalists have succeeded in Bosc- 
condensing atomic 52 Cr with a magnetic dipole moment 
of d = 6/ib [SI H], atomic 164 Dy with a magnetic dipole 



moment of d = 10/ie [IS] and atomic 168 Er with a mag- 
netic dipole moment of d = 7^b [16] . where /j,b is the 
Bohr magneton. Additionally, progress is being made, for 
example, towards the condensation molecular RbCs |17j . 
where the dipolar effects should be considerably larger 
than in the atomic 52 Cr, 164 Dy and 168 Er BECs. Clear 
dipolar effects have been observed in the 52 Cr, 164 Dy and 
168 Er BECs, though, in spite of relatively weak dipole- 
dipole interactions jTSJ IH] • 

For polarized dipoles, the anisotropic nature of the in- 
teraction leads to attraction when the dipoles are aligned 
head-to-tail (in the direction of polarization) and thus to 
an energetic instability in the homogeneous dipolar BEC 
(DBEC) [2U]. However, confinement in the direction of 
polarization can significantly stabilize the dipolar system 
against instabilities due to inelastic collisions [5T] and 
three-body loss processes [22]. To this end, the quasi- two- 
dimensional (q2D) geometry is sought after to suppress 
the attractive part of the ddi. Interestingly, a roton- 
maxon character is predicted to exist in the quasiparticle 
dispersion of the single-component DBEC in this q2D ge- 
ometry [2"5H2"5] . similar to that in superfluid 4 He, though 
from a different microscopic origin. In the q2D dipolar 
BEC, the roton-maxon dispersion can exist in the dilute, 
uncorrelated state, while it is precisely the correlations 
that gives rise the the roton in the superfluid 4 He sys- 
tem [2"o] . 

In this article, we consider a two-component BEC with 
both short-range and dipole-dipole interactions in the 
q2D geometry. Such a system can be realized by dual- 
condensing different atomic or molecular species, or by 
using different magnetic sublevels of the same atomic 
species [37] . Indeed, we find that roton physics manifests 
itself in a unique way in this two-component dipolar sys- 
tem, resulting in a first-order phase transition from the 
miscible to the immiscible state due to unstable roton-like 
quasiparticle fluctuations of the two-component BEC. In 
practice, this corresponds to a critical length scale at the 
IMT threshold that is set by the roton wavelength, be- 
ing on the order of the length scale of the tight trapping 
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potential. This is in stark contrast to the threshold tran- 
sition length scale of the non-dipolar condensate, being 
the phonon length scale, which is typically the longest 
available length scale in the system. 

Whereas other studies have characterized the IMT 
of dipolar Bose gases via full mean-field simulations in 
one dimension |28j and Monte-Carlo simulations in two- 
dimensions |29| , we perform full mean- field simulations in 
addition to an analytic Bogoliubov treatment, allowing 
us to characterize the IMT very efficiently in a large pa- 
rameter space. Additionally, we note that other theoret- 
ical works on non-dipolar binary BECs have predicted fi- 
nite wavelength phenomenon regarding, for example, the 
boundaries of immiscible systems [30H33] and quenches 
deep into immiscible parameter space [34] . These phe- 
nomena, however, are not roton-like in that the thresh- 
old defining the transition parameters is determined by 
long-wavelength, or phonon-like excitations in each case. 
The q2D binary DBEC is unique in that the IMT has 
finite wavelength character, or roton character, at the 
threshold, which is tunable as a function of many system 
parameters, including the tilt of the polarization field. 

This paper is organized as follows. In section [ill B| wc 
introduce the mean-field formalism for the interacting 
two-component BEC, discussing the three-dimensional 
homogenous case in section [lI A| and the quasi-two dimen- 
sional homogeneous case in section [Tl B| In scction [Tl C| we 
derive the two-component quasiparticle dispersion, which 
we use to identify both the familiar long-wavelength and 
the roton IMT in the large parameter space of the sys- 
tem. We present our results in section |III| including 
explicit solutions to the coupled Gross-Pitaevskii equa- 
tions, describing the condensate wave functions of the 
two-component BEC, for free and radially trapped ge- 
ometries. We conclude in section HVl 



II. MEAN-FIELD THEORY 

We consider an ultracold, dilute two-component Bose 
gas in the presence of two-body s-wave and dipole-dipole 
interactions, labeling the components by the indices j 
and k. The contact interactions between components j 
and k are characterized by the s-wave scattering lengths 
aj k through the pscudopotential 

V^{v-v')=g jk 8{v-v') (1) 

where gj k — 2ith 2 a.j k /Mj k is the contact interaction cou- 
pling matrix, Mj k — MjM k /(Mj + M k ) is the reduced 
mass matrix of two-body system and Mj denotes the 
mass of a boson of component j. Note that the in- 
traspecies (diagonal) reduced mass is just Mjj = Mj/2. 
Here, an = 021 is the interspecies s-wave scattering 
length and an and 022 are the intraspecies s-wave scat- 
tering lengths. The dipole-dipole interaction (ddi) po- 
tential for polarized dipoles is given by the potential 



Vf k (r) 
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where dj is the magnitude of the dipole moment of com- 
ponent j, d is the direction of the dipole polarization (as- 
sumed the same for both components) and r is the vector 
separating the two dipoles. The full two-body interaction 
potential is given by V}fc(r) = V? k {v) + V^ k {v). We can 
write the energy functional for the fully condensed sys- 
tem in terms of the condensate order parameters, or wave 
functions, ^j(r), 



E [{¥,•}] = / dr 



J2 **( r )Mr)*;,(r) 
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Here, hj(r) is the single particle Hamiltonian 



^•(r) = - 2 ^V 2 + C/,(r) 1 



(4) 



and Uj(r) is the trapping potential of component j. We 
normalize each ^j(r) to / dr|\&j(r)| 2 = Nj where Nj is 
the number of particles in component j and the total 
number of particles is given by N — Ylj=i 2 fy. In the 
energy functional ([3]), the sum over j, k — 1,2 implies 
summing over all four combinations of these indices tak- 
ing on these values. The factor of 1/2 in front of this sum 
takes care of the double counting in the intraspecies inter- 
action terms j — k and of the double counting produced 
by summing both {j, k} — {1, 2} and {j, k} = {2, 1} for 
the interspecies interactions. 

The coupled of Gross-Pitaevskii equations (GPEs) for 
this two-component system are derived by requiring that 
small variations of E 1 ^^-}] with respect to tyj vanish, 



giving 



2M. 



-V 2 + [/,(!■)+ ]T g jk \V k (r,t)f 



k=l,2 



J2 J dr'Vf k (r-r')\y k (r',t)f 



k=l,2 



*i(r,t).(5) 



(2) 



Here, we have generalized to the time-dependent form 
of these equations. The time-independent forms are re- 
covered by asserting the time dependence %(r, t) = 
^j(r)e~ ltljt where fj,j is the chemical potential of com- 
ponent j. Solutions ^j(r), corresponding to the fully 
condensed ground state, are found by minimization of 
the corresponding energy functional ([3| , and the dynam- 
ics can be studied by direct numerical integration of (|5j. 
In practice, as we explain in more detail later, we mini- 
mize ([3]) by evolving ^ in imaginary time |35) . 
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A. Homogenous Three-Dimensional System 

In the homogeneous three-dimensional (3D) geometry, 
there is no trapping potential and the miscible ground 
state of the two-component system can be described sim- 
ply by the condensate order parameters \&j (x) = y^ 13 , 

where nf is the 3D condensate number density of com- 
ponent j. For the non-dipolar case, the IMT of the homo- 
geneous system is characterized simply by the parameter 
A, introduced in section [I] This characterization can be 
made more rigorously by a linear stability study of the 
miscible state. Such a study is performed within the Bo- 
goliubov approximation, which results in a quasiparticle 
description of the two-component BEC. In this case, a 
pair of quasiparticles describes in-phase and out-of-phase 
two-component modes. The out-of-phase modes play an 
important role in characterizing the stability of the mis- 
cible state. We save a detailed discussion of such a theory 
for the case of dipolar interactions in the q2D geometry 
for the following section |II B| where the physics of the 
IMT is more complicated than for the non-dipolar case 
or homogeneous 3D cases. 

When one or both of the components in the homo- 
geneous 3D geometry possesses a non- negligible dipole 
moment, the polarization breaks the angular symme- 
try of the mean-field interaction potentials. In the one- 
component DBEC, this results in a quasiparticle disper- 
sion that is phonon-like (linear in the long-wavelength 
limit) but anisotropic. In this case, the speed of sound 
c depends on the angle q between the wave vector q 
and the dipole moment, c — > c(# q ). Similarly, the disper- 
sions of the two-component DBEC depend on (9 q and are 
characterized by a A parameter that depends on # q |20j , 



A-(0 q ) 



Gll(# q )G22(#q) 



1, 



(G) 



which determines the IMT threshold. If A(6» q ) < for 
any angle q , the components of the 3D two-component 
DBEC are immiscible. In the above equation, the 
Gjk{0q) matrices characterize both the short-range s- 
wave interactions and the ddi between components j and 



G jfe (0 q )= 5jfc {l + e^[3cos 2 q -l]} 



(7) 



and ejf characterizes the strength of the ddi between 
components j and k, e^f = bjk/ajk, where bjk — 
2Mjkdjdk/3fi 2 is the dipole length characterizing the 
ddi between components j and fc. With Mi = M2, 
the interspecies dipole length (j ^ k) becomes the ge- 
ometric mean of the intraspecies dipole lengths, 621 = 
b\2 = V^n^22, an d is thus uniquely determined by the 
intraspecies dipole lengths. These dipole lengths are de- 
fined so that a single component polarized DBEC in a ho- 
mogeneous 3D geometry requires a positive s-wave scat- 
tering length a,jj > bjj in order to energetically stabilize 
the system. 



It is interesting to consider the case where one com- 
ponent possesses a dipole moment and the other is non- 
dipolar, say, 622 = 0. In this case, A(# q ) is always maxi- 
mized when 9q = tt/2, corresponding to two-component 
quasiparticle propagation in the direction perpendicular 
to the dipolar polarization. In this system, the IMT 
threshold, A(6* q ) = 0, is first crossed in the direction per- 
pendicular to the dipole polarization, # q = tt/2. With the 
exception of the angular dependence, the stability con- 
dition resembles that of the non-dipolar two-component 
BEC system. 



B. Homogeneous Quasi- Two Dimensional System 

Recently, much effort has been directed towards realiz- 
ing trapping geometries with very strong confinement in 
one direction. Indeed, such geometries significantly sta- 
bilize dipolar gases that are polarized in the direction of 
the strong confinement by suppressing the attractive part 
of their interactions. When the characteristic interaction 
lengths of the trapped atoms or molecules are much larger 
than the harmonic oscillator length, this system is effec- 
tively two-dimensional (2D). However, when the interac- 
tion lengths are sufficiently smaller than the harmonic os- 
cillator length, the system develops q2D character where 
the zero-point motion in the trapped direction is impor- 
tant in characterizing the interactions [2H Hi] and pure 
condensation occurs at finite, as opposed to zero temper- 
ature [36 . A possible realization of this geometry uses 
a retro-reflected laser to create a one-dimensional optical 
lattice potential [57] . 

We model the two-component Bose system in the q2D 
geometry with the trapping potential potential Uj(z) = 
\Mj0j 2 (z 2 + p 2 /X 2 ) where A = lo z /uj p > 1 is the trap 
aspect ratio. When A 3> 1 and the interactions are rela- 
tively weak compared to the trapping energy huj z , we use 
the single mode approximation (SMA) where the conden- 
sate wave functions are assumed to have the separable 
form 



%(x) = <Pj(p)Xj(z) 



(8) 



where <Pj(p) is the in-plane wave function normalized to 
Nj, J dp\ifj(p)\ 2 = Nj, and Xj( z ) is the axial wave func- 
tion. We take Xj( z ) to be a Gaussian normalized to unity 
with width lj = ^h/MjUj z , 



X 3 (z) 



exp 



21) 



(9) 



We assume that both components are trapped in a har- 
monic trap with the same frequency Li z , but allow for 
different masses so the axial wave functions Xj( z ) can 
have different widths. 

Even in highly oblate traps, the separable ansatz (JsJ) is 
not exact, except in the case of a non-interacting system. 
The interactions, even when weak, also drive the axial 
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wave functions Xj( z ) away from the Gaussian form. How- 
ever, the SMA that we use here significantly simplifies the 
problem at hand, allowing us to explore a larger region in 
parameter space, and captures the relevant physics of the 
system [35] • The SMA is particularly beneficial in that it 
allows us to reduce the problem to a set of 2D equations 
by analytically integrating out the z-dependence in the 
coupled set of GPEs. For the dipolar mean field terms, 
this amounts to calculating an effective q2D interaction 
potential, given by 



p') = dz dz'xUz')Vf k (T-r')x 2 Az). (10) 



We handle this expression by transforming into momen- 
tum space, and thus need the Fourier transforms of Vf k (r) 
and Xk( z )- Without loss of generality, we consider a po- 
larizing field d = cosaz + sinax, which describes dipoles 
that are all tilted by an angle a off of the z-axis into the 
x-direction. The momentum space interaction for dipoles 
with this configuration is given by |20) 



ird ti \ 4-7T J (k z cos a + k x sin a) 



k 2 



1 



and the transform of the axial density is given by 

f 



Fid [Xj(z)] = exp 



(11) 



(12) 



where .Fid is the ID Fourier transform operator. Thus, 
the effective q2D momentum space ddi is given by [25|l39] 



v%(K) = D 3k F 



kpljk 



(13) 



where l jh = ^ + Z|)/2, D jk = ^g-^ is the ddi 

coupling matrix and F(q) = cos 2 aFj_(q)+sm 2 (a)F^ (q), 
where 



Fj_(q) = 2 - Z^qe q erfc(g) 



and 



F||(q) = -1 



3^^e g2 erfc( ? ), 



(14) 



(15) 



and erfc(g) is the complimentary error function of q |39) . 
This result ( 15 ) can be generalized to describe a polariza- 



tion field that is rotated by an angle rj off of the a:-axis, 
by taking q x -> q d , where q d = \J q x cos 2 77 + q y sin 2 rj is 
the wave number in the direction of the polarization tilt. 

Similar calculations can be carried out for the kinetic, 
potential and contact interaction terms to yield a cou- 
pled set of GPEs that govern the in-plane wave functions 
fj{p) of the q2D system, 



ihd t tpj(p,t) 



--— V 2 
2Afj 9 



2A 2 



Mjuj 2 z p 2 



fe=l,2 



9jk 



i^l 



MpA)\ 2 + / d P 'vUp-p')\<p k { P ',ttf 



tpj(p,t). 



(16) 



We calculate the dipolar interaction in the last term in 



Eqs. ( 16 1 by employing the convolution theorem, 



/ dp'vf k (p- 



p')\ Vk {p',t)\ 2 =T^ [v%(k p )h k (k p ,t)} , 



where Tin is the 2D Fourier transf orm operator and 
n k (k p ,t) = J-2D [\<Pk(p, t)\ 2 ] ■ Eqs. (16) fully describe 
the two-component DBEC in the q2D geometry (in the 
mean- field framework), where the components can have 
different masses and interaction character. It is interest- 
ing to note that the off-diagonal elements of the ddi cou- 
pling matrix Dj k are determined uniquely by its diagonal 
elements while the off-diagonal elements of the s-wave 
contact interaction coupling matrix gj k are, in principle, 
independent of the diagonal elements. Physically, the 
latter are determined by the microscopic structure of the 
components and are tunable via magnetic Fano-Feshbach 
resonance [101 EI] • 



C. Linear Stability: Bogoliubov Theory 

An instructive case to consider is the pure q2D case, 
or the case with no radial trapping potential (A — » 00) so 
the system is homogeneous in the x-y plane. In the mis- 
cible state away from any instabilities, the condensate 



(17) wave functions can be written as <pj(p,t) 



where rij is the integrated 2D density of the j th compo- 
nent. In the immiscible state, however, the continuous 
translational symmetry of the system is broken. We can 
study this immiscibility-miscibility transition (IMT) by 
considering small deviations from the ground-condensed 
miscible states in the form of Bogoliubov quasiparticles, 
ifij(p,t) — > \pri. (1 + 6ipj(p,t)) e~ v J 4 , where S <?C 1 and 



ipj(p,t) = uje' 



(18) 



where Uj and Vj are the Bogoliubov particle and hole 
amplitudes, respectively, and obey the normalization 
l u j| 2 — \ v j\ 2 = 1 03- We derive a set of equations for the 
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frequencies uj, the Bogoliubov de Gennes (BdG) equa- where hjk = 
tions, by linearizing the Bogoliubov ansatz about the momentum, 
small parameter 6 in the GPEs (16). The same equa- 



tions can be derived from a second quantized theory, by 
diagonalizing the full two-component Hamiltonian in the 
Bogoliubov approximation |43j . These equations can be 
written in matrix form as 1441 



/ijfe(q) is a function of the quasiparticle 



%\i = UJU 



(19) 



where u T = [u\, Vi, U2, ua] and the two-component BdG 
Hamiltonian % is given by the 4x4 matrix 



hjk(q) = n k 



9 3 k 



2ttZ 



jk 



D jk F 



q/ 



jk 



V2 



(22) 



H = 



Bii Bu 
B21 B22 



(20) 



and the Bjk are 2x2 sub-matrices, and are given by 



2 A/, 



Sjf. + hjk hj k 

h 2 q 2 c _ t 
2M, °J fe '* 



(21) 



jk 



jk 



An algebraic diagonalization of the BdG equations ( 19 ) 
yields the two-branch dispersion of the miscible q2D sys- 
tem, 



4(q) 



Wl(q)+w|(q)±W(w?(q) 



;(q)) 2 + 



Aq 4 



M X M 2 



(23) 



where uij (q) are the single component Bogoliubov disper- 
sions 



(24) 



The two-component dispersion ( 23 1 is identical to that 



of the homogeneous system with contact interactions [43. 
but with full momentum-dependent interaction coupling. 
We solve for the BdG eigenvectors u numerically. From 
these eigenvector solutions, we identify the upper (+) and 
lower (— ) branches of the two-component dispersion (23) 
as corresponding to in-phase and out-of-phase modes, re- 
spectively. As is clear from the form of the dispersion, 
u;_ will always be lower in energy than u + and is thus 
the relevant branch regarding stability of the miscible 
state. In this work, we consider only positive intraspecies 
and interspecies s-wave scattering lengths. In this case, 
and in the absence of the ddi, the single components 
are always stable and any instability, corresponding to 
Im[w_] 7^ 0, signifies a transition to an immiscible state. 
When the ddi is present, however, the single particle dis- 
persions can present dynamical roton instabilities when 
the effective dipole length is sufficiently larger than the 
positive s-wave scattering length, due to the momentum- 
dependence of the ddi in the q2D geometry [231 • As 
a result, imaginary parts of the lower branch of the two- 
component dispersion could correspond to transitions to 
immiscible states or to collapse of the miscible system, 
much like the roton collapse seen in the single-component 



DBEC. 

Before proceeding, we consider the most radical devi- 
ation from the simple Gaussian z-dependence expressed 
by Eq. ([9]), being phase separation in the trap (z) direc- 
tion. In the absence of the trap, miscibihty in this (z) 
direction corresponds to stable two-component quasipar- 
ticle propagation in this direction, which corresponds to 
stable quasiparticle propagation at an angle a off of the 
z-axis, where a is the polarization angle, for the homo- 



geneous 3D case discussed in section II A Thus, for the 
system to be miscible in the trap direction, the parame- 
ters must satisfy (from Eq. (§) 



A(a) > 0. 



(25) 



Indeed, all of the examples we consider here satisfy this 
criterion, so the assumption that there is no spatial sep- 
aration or immiscibility in the z direction is a good one. 
Even if Eq. (25) is not satisfied, however, the system is 



likely still miscible in the trap direction due to the ten- 
dency of the strong trapping potential to force a system 
to miscibihty. In this case, the criteria for miscibihty in 
the trap direction in the homogeneous q2D geometry fol- 
lows directly from previous studies of the IMT in a two- 
component BEC in a trap, as the ddi can be treated by 
using a modified s-wave scattering length in this case [24] . 
This was recently considered for specific dipolar species 
in ESI. 
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III. RESULTS 

For simplicity, we proceed by considering a two- 
component system with equal integrated densities n\ — 
n 2 = n and masses Mi = M 2 = M, which we refer 
to as the "balanced" system. Having such a balanced 
system results in the two components sharing an axial 
wave function, Xi( z ) = X2( z )- In the homogeneous 3D 
geometry, the stability of the mixture does not depend 
on the direction of the dipole moment. Such is not the 
case, however, in the q2D geometry. If the dipoles are 
polarized in-plane, corresponding to a = tt/2, the sys- 
tem acquires the same energetic stability criterion as the 
homogeneous 3D case. If the dipoles are polarized in 
the trap direction, corresponding to a = 0, the criterion 
for energetic instability becomes a,jj > —2bjj [53]. For 
now, we consider the a = case. Additionally, from here 
forward we set nl^ = 1, so that the interactions are char- 
acterized solely by the ratio of the interaction lengths to 
the axial harmonic oscillator length l z . 



A. Roton Immiscibility for a = 

For the short-range s-wave interactions, we consider 
scattering lengths that would result in a slightly immis- 
cible system in the absence of the ddi. For now, we take 
dli/h — CL22/I7, — 1-0 and ayijlz = 1.05. Additionally, 
we fix h\\jl z = 2/3 and explore the stability of the mis- 
cible system as a function of 622/^ using the Bogoliubov 
theory laid out in the previous section. 

It is straightforward to predict the behavior of the two- 
component system when bn = 622, that is, when both 
components are equally dipolar. In this case, the intra- 
and inter-species dipolar interactions are equally repul- 
sive and their effects cancel in the dispersion cj_(q). As 
a consequence, the ddi plays no role in the IMT of the 
system. The IMT is instead solely determined by the 
short-range parameters. Indeed, this is seen in figure [l] 
where the imaginary part of the lower branch of the two- 
component Bogoliubov dispersion, w_(q) from Eq. ( |23[ ), 
is shown £iS 0jS 0/ function of the wave number ql z . For a 
range of &22A2 — bu/l z = 2/3, there is a long-wavelength 
immiscibility, labeled and shown by the shaded region in 
this figure. For 6 2 2 7^ ^11 , however, the difference in the 
dipole moments of the two components plays an impor- 
tant role in characterizing the stability of the miscible 
state. For example, reducing 622 relative to b\\ stabilizes 
the long-wavelength instability of the miscible state, as 
seen in the stable gap for a range of 622 /h m figure [T] 

Another important example of this, and indeed a key 
result in this paper, occurs for smaller values of 622- For 
the parameters given above, and for 622A2 ^ 0.1, there 
exists another region where Im[w_(q)] ^ 0. This region 
occurs at finite, nonzero wave number, and is charac- 
terized by the softening of a roton-like feature in the 
cj_ quasiparticle dispersion. As in the case of a single 
component DBEC, the softening of the roton dispersion 
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FIG. 1: (color online) The phase diagram of the bal- 
anced two-component system in the q2D geometry with dipole 
lengths b\\jl z — 2/3 and scattering lengths an = a22 = 1 and 
ai2 = 1.05 as a function of the dipole length of component-2, 
622/^2 and quasiparticle wave number ql z . The familiar long- 
wavelength immiscible region and the roton-immiscible region 
are labeled and colored, and the white region signifies the ex- 
istence of a stable miscible state. The roton-immiscible region 
extends down to 622/^z = 0. The contours in the shaded re- 
gion mark increments of 0.1lj z . 



signals an instability of the q2D-homogeneous system. 
Unlike the single component DBEC, however, the insta- 
bility does not lead to collapse but instead results in an 
immiscible density pattern. The features in the density 
pattern have a length scale that is the inverse of the mo- 
mentum at which the the roton dispersion touches the 
axis. We term this kind of phase-separation 'roton im- 
miscibility' While the case b 2 ifl z = is not shown in 
this figure (due to the logarithmic scaling), the region of 
roton immiscibility extends down to this limit. 

The limit of = is, in fact, useful in reveal- 

ing the nature of the roton immiscibility. In this case, 
the only dipolar interactions in the system occur within 
component 1. The dispersion of component 1, oj 1 (q), 
does not possess a roton minimum, but is not purely 
phonon-like, either, due to the momentum dependence 
of the interaction /in(q), Eq. ( p2| . This dispersion is 
shown by the dash-dotted line in figure [2] along with the 
dispersion of component 2, cj2(q) (dashed line), which is 
purely phonon-like due to the purely short-range nature 
of these interactions. The a\ijl z ratio that characterizes 
the interspecies repulsion is then sufficiently strong to 
drive the system to immiscibility, but on a length scale 
set by the emerging roton in component 1. The lower 
branch of the two-component dispersion, w_ (q) , is shown 
in figure [2j as well, where its imaginary part, signifying 
the transition to immiscibility, is labeled and shaded. 

We map the roton immiscibility phase boundary for 
the balanced system with an/l z — 022/^2 = 1 and 
b 22 /l z = in figure [j^a). This figure shows the IMT 
threshold as a function of the interspecies scattering 
length a\2jl z and the dipole length of component 1, 
b\i/l z - Interestingly, the roton immiscibility persists for 
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FIG. 2: (color online) Dispersions of the balanced, misci- 
ble, two-component q2D DBEC with an/l z = ai%/l z = 1, 
a\ijl z = 1.05, bn/l z — 2/3 and b2i/l z = 0. The lower branch 
of the two-component dispersion, w_(q), corresponding to 
out-of-phase quasiparticle excitations, exhibits a roton-like in- 
stability corresponding to a transition to immiscibility. The 
imaginary part of the dispersion is shaded in teal. 



all non-zero values of b\\ when the intraspecies scatter- 
ing length is larger than a critical interspecies scattering 
length. This is shown in figure [3]^a), where we plot the 
threshold intraspecies scattering length as a function of 
bii/l z , above which the system is immiscible (shaded) 
and below which the system is miscible. Interestingly, for 
^ 0.6 the threshold for roton immiscibility occurs 
at a\ijl z < 1. In the absence of any ddi, the transition 
to miscibility occurs when a\ijl z > 1, which is shown by 
the black dotted line in this figure. 

The onset of instability at a finite momentum is a 
roton-specific feature that signals a first-order (zero- 
temperature) phase transition. The inverse of the critical 
momentum g cr jt at which the roton in the w_ disper- 
sion softens indicates the length scale at which immis- 
cible (single condensate) density features nucleate if the 
system evolves from a homogeneous ground state. Fig- 
ure [3jb) shows the critical wave number q c ^itl z at which 
the dispersion w_ (q) first develops a non-zero imaginary 
part with increasing a^/h as a function of bn/l z . For 
bii/l z = 0, corresponding to a completely non-dipolar 
system, the immiscibility transition occurs at q C ritl z = 0, 
which is the familiar long-wavelength immiscibility found 
in BECs with only short-range interactions. As b\\jl z is 
increased, however, the critical wave number for immisci- 
bility increases, corresponding to a transition to immisci- 
bility with roton, or finite wavelength character. Indeed, 
when bn/l z = 1, the transition to immiscibility occurs 
at the large wave number ql z ~ 2, corresponding to a 
transition wavelength A cr it ~ ttI z . 

We proceed to further identify the roton instability of 
the miscible state as a transition to immiscibility, and not 
an instability to a collapsed state, by directly solving the 
coupled GPEs, Eqs. (16), in the absence of a radial trap, 



FIG. 3: (color online) Critical interspecies s-wave scattering 
length a^/lz (top panel) and wave number (bottom panel) for 
the balanced two-component q2D DBEC as a function of the 
dipole length of component 1, bn/l z , with an /l z = 0.22/1? = 1 
and &22/L = 0. The parameters for which the system is 
immiscible are shaded teal in the top panel. The non-dipolar 
IMT is seen for bu/l z = 0, where the IMT threshold occurs at 
aw/lz = 1 (A = 0) and q C xnl z = 0. The roton immiscibility 
emerges for nonzero 6n/i z , and decreases in characteristic 
wavelength as b\\jl z is increased. The black "x" in the top 
panel marks bn/l z = 2/3 and a\2fl z = 1.05, the parameters 
for the example roton IMT used in figures [2] and [5] 



state, we sample the condensate wave functions of the two 
components on a numeric grid and employ the imaginary- 
time evolution algorithm, stopping when the total energy 
is converged to a part in 10 8 . We choose a grid of size 
512 x 512 with the spatial extent x G [— x max , x max ] and 
V G [-J/max, 2/max]- These real-space limits x max and y max 
are chosen so that the spatial resolution of the numeric 
grid is sufficiently smaller than both the healing length 
of the system and the axial harmonic oscillator length 
l z . To initiate the algorithm, we break the symmetry of 
the system by seeding the initial homogeneous guesses 
for the condensate wave functions with numeric noise in 
the form of two-component quasiparticles. That is, we 
take for component j [46] 



cpj(p,t = 0) -> v^jl + ^v^e" 

x [uj.ie-^P + vl l e l ^ p ]} (26) 

where {a{\ are random numbers such that cv; G [0,1] 
for all I and V\ = Ni/Nq, Nq — nl 2 z and Ni is the number 
occupation of the quasiparticle state of component j with 
energy given by the Bose-Einstein distribution, 



1 



r. 



(27) 



where T is the temperature of the system. While we 
choose T = 100 nK, this does not carry strong physical 



corresponding to A — > 00. To find the stationary ground meaning when evolving Eqs. (16) in imaginary time, as 



this equation is dissipative. We calculate the Ujj and 
v* i quasiparticle amplitudes via numeric diagonalization 
of (21). Additionally, we choose q; such that an integer 



Component 1 



Component 2 



number of quasiparticle wavelengths fit in the grid, so 
that the periodic boundary conditions of the system are 
satisfied. Such a restriction, however, is not important 
when evolving the system in imaginary time, as any un- 
physical, high-momentum components of the initial wave 
functions that are rooted in relaxing this restriction will 
quickly dissipate. Indeed, we choose a random direction 
for each q; to introduce noise that does not share the 
symmetry of the numeric grid. Finally, we sample uj-.i 
from the dispersion w_(qz), noting that these energies 
are the same for both components. 



Using this algorithm to find stationary solutions of the 
coupled GPEs, Eqs. ( |16[ ), we first explore the case when 
bii/l z = = 2/3. As mentioned above, this case is 

immiscible but the dispersion is phonon-like and does not 
result from the softening of the roton dispersion feature. 
In accordance, the ground state solution to the coupled 
GPEs, shown in figure |1J exhibits immiscibility at the 
longest available wavelength. 



The case where &22A2 = and bn/l z = 2/3, however, 
is in stark contrast to the case where both components 
possess equal dipole moments. The two-component Bo- 
goliubov theory predicts a roton instability of the misci- 
ble mixture in the w_(q) dispersion. Indeed, we find that 
the ground state of this system exhibits immiscibility on 
a much shorter length scale, being the length scale of the 
roton. The condensate densities of the two components 
are shown for this case in figure [5| 

For the cases shown in figures [3] and [5] the ground state 
solutions depend on the numeric noise, in this case the oc- 
cupation of single quasiparticles with well-defined, albeit 
random, standing wave orientations, that is used to seed 
the condensate at the beginning of the imaginary time 
evolution. For example, if wave vectors are chosen such 
that q = qx, we find that the immiscible ground states 
possess density fluctuations only in the x-direction. Sim- 
ilarly, we performed simulations with random noise in 
the form of small-amplitude random numbers sampled 
at each grid point. The immiscible ground state for the 
case with one dipolar and one non-dipolar component 
(the roton immiscible ground state) shows a speckle, or 
bubble-like pattern, fluctuating on a length scale ~ Al z . 
In the cases discussed here, the origin of the sensitivity 
to the symmetry of the initial seeding of the wave func- 
tions is the fact that there is no internal system bias for 
the direction of the immiscibility (a = 0). However, the 
anisotropy of the ddi can be exploited to introduce an 
anisotropic momentum-dependence in the system inter- 
actions (a ^ 0), thus breaking the azimuthal symmetry 
of the system and creating a directional bias. 





FIG. 4: (color online) Densities n 1 (p) = \ipi(p)\ 2 (left 
panel) and «2(p) = \(p2(p)\ 2 (right panel) corresponding to 
stationary solutions of the coupled GPEs ( 16 1 (converged 
in energy to a part in 10 s in imaginary time evolution) for 
the non-dipolar balanced system, with au/l g = (122/lz = 1, 
o-u/lz = 1-05 and bn/l z — 622/^ = 0. Here, the immisci- 
bility is phonon-like and occurs on the longest length scale 
available in the system. 
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FIG. 5: (color online) Densities rii(p) = \(pi(p)\ 2 (left panel) 
and Ti2(p) = \(f2(p)\ 2 (right panel) corresponding to station- 
ary solutions of the coupled GPEs (161 (converged in energy 



to a part in 10 in imaginary time evolution) for the balanced 
system, with an/l z = 022/^2 = 1, a^/lz = 1-05, 622/^ = 
and bii/l z — 2/3. Here, the immiscibility is roton-like and 
occurs on the length scale ~ 2irl z . 



B. Roton Immiscibility for a > 

The behavior of the quasiparticle dispersion of the 
single-component q2D DBEC was found to exhibit in- 
teresting quasiparticle dispersion character as a func- 
tion of polarization direction . More specifically, an 
anisotropic roton emerges as a function of polarization 
angle a for certain interaction strengths and densities, at 
wave vectors perpendicular to the direction of the polar- 
ization tilt. Here, we find an analogous phenomenon in 
the two-component system. 

For now, consider a two-component BEC with the pa- 
rameters used above, but with a smaller dipole length in 
component 1, bn/l z = 1/4, and no dipolar character in 
component 2, 622/^z — 0- If the polarization axis is per- 
pendicular to the trapping direction, corresponding to 
a = 0, the ground state of this system is a miscible mix- 
ture of the two components. As the polarization axis is 
tilted, however, an instability to a miscible state emerges 
as the character of the ddi in component 1 becomes more 
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FIG. 6: (color online) The shaded region shows the non- 
zero imaginary part of the two-component dispersion oj_ (q) 
in the direction perpendicular to the polarization tilt for 
aii/h — a-22/lz — 1, au/h = 1-05, bn/l z — 1/4 and 
biijlz = as a function of polarization angle a. and quasi- 
particle wave number ql z . The miscibility of the system for 
a = is reflected in the purely real character of the dispersion 
at this polarization angle, and the roton-immiscibility is seen 
as the emergence of a non-zero imaginary part of t<j_(q) at 
finite wave number ql z ~ 1 at a/ir ~ 0.05. The contours in 
the shaded region mark increments of 0.1w z . 
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FIG. 7: (color online) The critical polarization angle a (top 
panel) and the critical wave number q CI itlz (bottom panel) 
for the balanced two-component q2D DBEC with an/l z = 
a-22jl z = 1 and 622/^2 = for various dipole lengths of com- 
ponent 1, bn/l z , as a function of the interspecies s-wave scat- 
tering length a^/lz- The immiscible parameters are shaded 
in the top panel. Here, the roton immiscibility emerges as a 
function of polarization angle for q cr i t perpendicular to the 
direction of tilt. 



attractive in this direction. We plot the imaginary part 
of the lower branch of the two-component dispersion rela- 
tion, uj_ (q) , for the component of the quasiparticle wave 
vector that is perpendicular to the polarization tilt as a 
function of the tilt angle a and the quasiparticle momen- 
tum in figure [6j The miscible ground state is reflected in 
the purely real dispersion at a = 0. A non-zero imagi- 
nary part emerges at a/-7r ~ 0.05 at finite, non-zero wave 
number, corresponding to a roton immiscibility that is 
emergent with polarization tilt. This roton immiscibil- 
ity has the same origin as that discussed earlier, being 
the momentum dependence of the ddi in the q2D geom- 
etry, though this momentum dependence is anisotropic 
and possesses angular dependence for a 7^ 0. Indeed, 
no immiscibilities are predicted in the direction of the 
polarization tilt for any a, suggesting that the striped 
structure of the immiscible state can be controlled by 
the proper adjustment of the polarization field. 

In figure [7[ we extend this result to characterize the 
onset of roton immiscibility as a function of tilt angle 
a and interspecies scattering lengths for various b\\/l z . 
In figure Qa) , the critical scattering length for the tran- 
sition to immiscibility is shown as a function of the in- 
terspecies scattering length for various dipole lengths of 
component 1, b\\jl z — 1/4,1/2,2/3. Again, we empha- 
size that component 2 is non-dipolar, so 622/^2 = here. 
In figure [7^b) , the corresponding critical wave number 
9criJz is shown. For all dipole lengths, the onset of im- 
miscibility as a function of ayijl z occurs at q cr itl z = for 
a critical polarization angle cv cr it = 7r/2, corresponding 
to the dipoles being polarized in the plane of symmetry. 
This is expected, though, as the interaction (22) exhibits 
no momentum-dependence for a = it/ 2 in the direction 



perpendicular to the polarization tilt and the criteria for 
miscibility can thus be extracted from that of the homo- 
geneous two-component DBEC, given in Eq. Q. Indeed, 
the values of a 12 /l z for the long- wavelength IMT thresh- 
olds at g cr it^z = can be found by setting A(ir/2) =0 and 
solving for ay2,jl z . As a,Yijl z is increased beyond these 
long-wavelength threshold values, however, the critical 
polarization angle decreases and eventually approaches 
«crit = 0, which characterizes the roton immiscibility dis- 
cussed earlier. In this parameter range, q cvl tl z > 0, and 
the critical wave number approaches characteristic roton 
wave numbers for larger a\2jl z . 

To further demonstrate the role that the polarization 
direction plays in the roton immiscibility of the two- 
component q2D DBEC, we model a time-dependent pro- 
cess (perhaps an experimental scenario) in which the po- 
larization field a is tilted as a function of time, driv- 
ing the system to an immiscible state through manip- 
ulating the anisotropy of the ddi in component 1. We 
model this scenario via direct time-dependent integration 
of the coupled GPEs, Eqs. (16). We take the param- 



eters for the balanced system introduced earlier, with 
a\x/l z = a 2 ijl z = 1, 0,12/h = 1.05, bn/l z = 1/4 and 
^22/^2 = 0. To begin, we consider a = at time tu> z = 0, 
corresponding to a miscible system. To break the sym- 
metry of this homogeneous ground state, we seed the 
condensates with quasiparticles as given in Eqs. (26) 
and (27), taking T = 100 nK. Then, we linearly ramp 



the polarization angle to a final value ahold = tt/4 over 
a time i ra mp, we hold the polarization angle at ahold for 
a time £hoid> a- n d we ramp the polarization angle back to 
a = over the time t ramp . For this simulation, we take 
^ramp^2 = 200 and ihoid^z = 400, so that the time scales 
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of the polarization tilt are much greater than the other 
time scales in the system, being the inverse quasiparticle 
energies, and the transition is thus to a good approxima- 
tion adiabatic. 

Because the real-time evolution of the coupled GPEs 
preserves total energy and is non-dissipative, we seed the 
condensate with quasiparticles with wave vectors in the 
x or y directions only. This ensures that the periodic 
boundary conditions of our numeric grid are satisfied and 
there are no unphysical high-momentum fluctuations at 
the edges of the system. As a result, the initial quasi- 
particles propagate only in the x and y directions. Thus, 
choosing a tilt direction of rj = or rj = tt/2 (where 77 
is the tilt angle relative to the x-axis) would share the 
symmetry of the initial quasiparticle noise. Instead, we 
choose T) — tt/5. According to the Bogoliubov theory 
presented above, the corresponding immiscibility should 
result in density stripes along the direction rj = tt/5. 

We present results from this simulation in figure [8j 
where we plot the densities of components 1 and 2 as a 
function of time. Column (a) shows the densities at time 
tuj z = 0, where the densities are approximately equal, 
aside from the initial seeding that is not visually notice- 
able in the given contour scaling. Column (b) shows the 
densities at time tu> z — 120, when the polarization angle 
just exceeds a C rit and the system begins to exhibit immis- 
cible character. Column (c) shows the densities at time 
tcj z — 200, at the end of the polarization angle ramp, and 
when the system exhibits full immiscibility. The immisci- 
bility clearly has roton character, and as anticipated, the 
direction of the immiscible stripes coincides with the tilt 
direction, at rj — n/5. Columns (d) and (e) show the den- 
sities at times tuj z — 760 and tuj z = 1200, respectively. 
These columns demonstrate the mechanism for symme- 
try breaking of the striped roton immiscible state, as the 
polarization angle tilts back to a = 0. While the system 
remains immiscible, the stripe character is completely 
lost by tu) z = 1200. Here, the introduction of strong 
phase fluctuations during the transition to immiscibility 
prevents the system from returning to the miscible state 
over the time scales considered here. 

The issue of thermalization in multicomponcnt BECs 
is of increasing interest with the finite temperature stud- 
ies of F = 1 spinor condensates, such as the 87 Rb sys- 
tem [57H1S]. Recently, it was shown that such systems 
also do not thermalize over even very long time scales, 
following a quench from the polar to the ferromagnetic 
state [501. 



C. Radially Trapped Case 

Until now, we have considered only the case where the 
q2D system is homogeneous in the x-y plane. In a re- 
alistic experimental scenario, the trapping potential will 
have a finite radial extent. For the roton immiscibility 
to persist in this geometry, we expect that the trapped 
system must be such that the density of the components 



is sufficiently large over a transverse length scale that is 
sufficiently greater than the roton length scale. To in- 



vestigate such a claim, we solve the coupled GPEs ( 16 ) 



numerically in the presence of a radial trap with aspect 
ratio A = 10, where A = uj z /oj p . While this trap aspect 
ratio is seemingly small, it serves to model a larger trap 
aspect ratio, as our solutions are free to extend in the 
radial direction but are Gaussians with a fixed width in 
the axial (z) direction. 

Much like the non-dipolar system, we find that a 
larger interspecies repulsion (corresponding to a smaller 
A) is necessary for the immiscibility to occur in a ra- 
dial trap [2]. Thus, we take a\\/l z = 022/^2 = 1 and 
a i2 /lz = 1-1 with, as before, 622/^z = 0, and investigate 
solutions for varying b\\jl z . We solve the coupled GPEs 
via imaginary time propagation, where we take the ini- 
tial wave functions to be Gaussians with small amplitude 
random noise sampled on the numeric grid. We converge 
the energy of each component to a part in 10 6 . 

We present some results in figure [9] for (a) bn/l z = 0, 
(b) bn/lz = 1/4, (c) bn/l z = 1/2 and (d) b n /l z = 2/3. 
For b\\jl z = 0, the system is non-dipolar and the IMT is 
characterized by the contact interaction strengths alone. 
In this case, the immiscibility is seen as the splitting of 
the components at the largest available length scale, be- 
ing the radial extent of the system. For b\\jl z = 1/4, 
the system is still immiscible, but the ddi is not suffi- 
ciently strong to induce roton immiscibility. Instead, the 
system exhibits a long wavelength immiscibility, but now 
the density of component 1 is pushed to the outside of 
the trap due to the intraspecies interactions of compo- 
nent 1 being greater than those of component 2. The 
roton immiscibility emerges near b\\jl z — 1/2, as we see 
in row (c) of figure [9j Because we initiate the imagi- 
nary time evolution with randomly sampled noise on the 
numeric grid, the system does not prefer a direction to 
break the symmetry of the miscible state, as is shown for 
the homogeneous q2D system in figure[5j Here, the roton 
immiscibility manifests as the dipolar component forms 
"bubbles" with spacings on the order of 10Z 2 , just larger 
than the characteristic roton wavelength A ro ton = 2ttI z . 
For b\\/l z = 2/3, shown in row (d), however, the mean 
spacing between clumps is on the order of A rot0 n- 

Because of the effectively stronger self-repulsion of 
component 1 (due to the ddi), it possesses a finite den- 
sity at radial extents greater than component 2. This 
is seen in all cases where b\\jl z > 0, in rows (b), (c) 
and (d) of figure [9] Interestingly, for smaller condensate 
densities, we find that increasing bn/l z can result in the 
majority of component 1 being pushed to the outside of 
the trap, suppressing the roton immiscibility. For inter- 
mediate b\\/l z , however, the roton immiscibility persists 
near the boundary of the two components. We find that 
such a phenomenon exists in more oblate traps, as well, 
where the radial trapping potential does not force the 
components to overlap in high density regions. 

It is interesting to note that the results in figure [9] are 
very similar to those presented in |27j . where a dipolar 
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FIG. 8: (color online) Condensate densities rii(p,t) = \<pi(p, t)\ 2 (top row) and ri2(p,t) = \(p2(p,t)\ 2 (bottom row) as a 
function of time, where the polarization angle is tilted from a value where the system is miscible (a = at t = 0, column (a)) 
to a value where the system is roton immiscible (ahold = 7r/4) at an angle 77 = 7r/5 off of the a>axis. In (a), the system is clearly 
miscible. Over a time t ramp oj z = 200, a is linearly ramped from a = to Qhoid = 7r/4. It is held at «hoid f° r thoid^z = 400, 



then is linearly ramped back to a — over t T 



200. 



In column (b), the densities are shown at time tuj z = 120, where 
the system is just starting to exhibit immiscible character. The densities are shown at time tui z = 200 in column (c), where 
the system is fully immiscible on the roton length scale. The densities are shown for toj z = 760 in column (d), as the angle has 
returned to a value corresponding to miscibility, and column (e) shows the densities at time tu z = 1200, where the polarization 
angle has returned to a — but they system is not rethermalized due to the introduction of strong phase fluctuations in the 
process of spatial seperation en route to miscibility. 



component and a non-dipolar component are separated 
in the direction of polarization by an external magnetic 
field gradient, forcing immiscibility in this direction. For 
a sampling of interaction strengths and field gradients, 
pattern formations on finite length scales, reminiscent 
of a classical magnetic ferrofluid, are predicted by di- 
rect solutions to the coupled GPEs that describe the sys- 
tem. These patterns are similar to those seen in figure |9j 
though perhaps from a different physical origin, and de- 
pend strongly on noise that is used to break the symme- 
try of the initial guess for the condensate wave functions. 

As is clear from our preceding discussion, the roton 
immiscibility depends strongly on the interaction param- 
eters in the system under consideration, and special care 
must be taken in proposing a candidate species with 
which to investigate such a phenomenon. One such can- 
didate is the ground 7 53 state of atomic 52 Cr prepared 
in the mj = spin projection for the non-dipolar com- 
ponent and the mj — —3 spin projection for the dipo- 
lar component, which has a magnetic dipolc moment 
of d = 6/ig where [ig is the Bohr magneton. While 
the relaxation lifetimes of the mj > —3 projections of 
the 7 63 state of 52 Cr were measured to be relatively 
long |51j . spin-exchange collisions, occurring on a time 
scale ~ .Is, may limit the experimental feasibility of us- 
ing 52 Cr. Nevertheless, for a trap with an axial frequency 
of lu z = 2ir x 2 kHz, the lifetime of the roton immisci- 
bility, for the maximum imaginary part of w_ (q) being 
~ 0.1u) z , is characteristically ~ 0.80 ms. The critical in- 
tegrated density for roton immiscibility in this case is 



n 2D 4.0 x 10 cm , corresponding to a maximum 
3D density of ~ 7.2 x 10 15 cm~ 3 . The immiscibility 
for the case seen in figure [9^c) with trap with axial fre- 
quency lo z — 2tt x 2 kHz would require a total number of 
52 Cr atoms on the order N ~ 6.2 x 10 6 with scattering 
lengths an = 022 = 30.4 an, and an = 33.4 ao, where an, 
is the Bohr radius. This is a case that is possibly achiev- 
able by sufficient manipulation of the Fano-Feshbach res- 
onances [41] . 

For more strongly interacting species, however, the ro- 
ton immiscibility emerges for a smaller critical density or 
particle number. Other possible dipolar species include 
atomic Dy, which has been recently Bose-condensed [TS] , 
and polar molecules, which can possess relatively large 
electric dipole moments. For the trap discussed above, 
the roton immiscibility emerges for a total number of 
164 Dy atoms of N ~ 400 x 10 3 , or a critical maximum 
density of ~ 2.6 x 10 15 cm~ 3 . An experimental study 
of the magnetic Fano-Feshbach resonances in atomic Dy, 
however, has yet to be performed. For polar molecules 
of, say, RbCs [T7] with an electric dipole moment of 
d = 0.5 Db, a critical density of ti 3 d ~ 9.0 x 10 13 cm~ 3 , 
or a critical particle number of N ~ 8900 molecules is 
needed. The realization of one dipolar and one non- 
dipolar component for the molecular case is unclear, how- 
ever, a mixture of, say, RbCs and Rb may be possible in 
the near future. Additionally, we have checked that ro- 
ton immiscibility exists for an appropriate set of s-wave 
scattering lengths in alkali atom (with ~ 1 fj,g magnetic 
dipole moments) and Cr or Dy mixtures. 
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IV. CONCLUSION 



The long-range and anisotropic nature of the ddi plays 
an interesting, nontrivial role in the physics of many- 
body systems. Here, we focus on the case of a two- 
component Bose-Einstein condensate and show that, for 
a set of specific interaction parameters and trap geome- 
tries, the system exhibits immiscibility with roton char- 
acter, where the immiscibility in a non-dipolar system 
or a homogeneous 3D dipolar system is strictly phonon- 
like. In particular, we find that the roton immiscibility 
occurs in the q2D geometry when the interspecies and 
intraspecies s-wave scattering lengths have comparable 
values, while one component is not (or negligibly) dipolar 
and the other component possesses a dipolar length that 
is comparable to the s-wave scattering lengths. By em- 
ploying the Bogoliubov theory to the q2D homogeneous 
two-component system, we calculated a two-component 
dispersion to efficiently characterize the stability of the 
miscible state and the parameters that define the IMT 
threshold for the familiar phonon-like and the roton im- 
miscibilities. Direct simulations of the coupled GPEs 
verify these results and reveal interesting dynamic and 
symmetry-breaking features of the immiscible phase. Re- 
cent experimental progress inspires confidence that the 
roton immiscibility may be observable in the near future. 



FIG. 9: (color online) Stationary densities ni(p) = |i^i(p)| 2 
(left column) and ri2(p) = \tp2{p)\ 2 (right column) of the 
balanced two-component DBEC in a trap with aspect ratio 
A = 10 for au/lz — ani/h = 1, oia/Z* = 1.1 and 622 //a = 0. 
The rows show solutions for various dipole lengths of compo- 
nent 1, being a) b lx /l z = 0, b) bu/l z = 1/4, c) bu/k = 1/2 
and d) b\\jl z = 2/3. In row (a), the system in non-dipolar 
and the interactions are thus balanced, resulting in an im- 
miscibility on the length scale of the trap. In row (b), the 
presence of the relatively small dipole moment in component 
1 breaks this symmetry, but the immiscibility is still long- 
wavelength and component 1 thus engulfs component 2 due 
to its greater self-repulsion. The roton immiscibility emerges 
for bn/l z > 1/2, as is seen in rows (c) and (d). While both 
immiscible states exhibit ordering near the roton length scale, 
the length scale is smaller in (d) where the dipole length is 
larger, as suggested by the results shown in figure [3] 
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